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Abstract. For chaotic systems there is a theory for the decay of the survival 
probability, and for the parametric dependence of the local density of states. This 
theory leads to the distinction between "perturbative" and "non-perturbative" 
regimes, and to the observation that semiclassical tools are useful in the latter 
case. We discuss what is "left" from this theory in the case of one-dimensional 
systems. We demonstrate that the remarkably accurate uniform semiclassical 
approximation captures the physics of all the different regimes, though it cannot 
take into account the effect of strong localization. 



1. Introduction 

The quantum mechanical state of a particle, or of a system, is represented by the 
probability matrix p. This object corresponds to a classical distribution p cl (Q,P) 
in phase space, where Q and P are the canonical coordinates of the system. One 
way to represent the probability matrix is by using the Wigner function p(Q,P). 
Many calculations in quantum mechanics reduce eventually to calculation of a trace 
over a pair of probability matrices. This includes in particular calculations of the 
local density of states (LDOS) JTj, and calculations of the survival probability V(t) 
0. Strongly related are calculations of Franck-Condon factors for non-adiabatic 
transitions between Born-Oppenheimer surfaces PJ- Lately El E] there is a 
special interest in calculation of "fidelity" (also known as "Loschmidt echo") in the 
context of quantum computation. This "fidelity" is in fact a synonym for a "survival 
probability" V(t) which is calculated for a particular dynamical scenario (namely, the 
forward evolution is followed by a reversed evolution with a perturbed Hamiltonian). 

The calculation of a trace over a pair of probability matrices has a clear classical 
limit. To be specific let us consider the survival probability, which is defined as 

V{t) = \{M^t)\ 2 = trace(ptpo)- (1) 

where pt and ipt are the probability matrix and the associated wavefunction of an 
evolving quantum mechanical pure state. In the Wigner representation the trace 
operation means dQdP / (2TrK) d integral over phase space, where d is the number 
of freedoms. It should be emphasized that the Wigner representation is quantum- 
mechanically exact. The classical limit /approximation is obtained by treating p(P, Q) 
as a classical distribution, whose evolution is governed by classical equations of motion. 

To take the classical limit as a leading order approximation for V(t) is one 
possibility. Another possibility is to use perturbation theory, in particular "Fermi 
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golden rule" (FGR). It should be obvious that the results that are obtained by using 
these methods are typically very different from each other. 

A major objective of recent studies jS] is to understand the limitations of 
perturbation theory on the one hand, and to explore the capabilities of the 
semiclassical tools, on the other. A specific question is whether the decay of V(i) 
is of perturbative nature (e.g. FGR type J), or else whether it is of classical nature 
("semiclassical" type). The main approach towards this question is to allow the 
specification of a control parameter that represents the "strength" of a perturbation. 
Depending on the value of of this control parameter one may have either perturbative 
or semiclassical behavior. The first publications that have taken this approach, 

led to the distinction between "perturbative" and "non-perturbative" regimes, and to 
the realization that the "semiclassical" behavior is contained in the" non-perturbative" 
regime. Later the idea was adopted into the context of "fidelity" studies. 

The above mentioned studies have mainly concentrated on quantized chaotic 
systems. The chaos assumption allows simplification of certain calculations. In 
particular one can invoke the random matrix theory (RMT) conjecture, in order to 
obtain some "generic" results. The natural question that arises is whether some of 
the general theory regarding the LDOS and V(t), applies also to the world of one 
dimensional (d — 1) systems. 

A central observation in the theory of quantized chaotic systems is the existence 
of two distinct energy scales. One is the mean level spacing (A oc h d ), while the other 
is the bandwidth (A& oc %). The latter is semiclassically related to the correlation time 
of the classical motion. The dimensionless bandwidth is defined as b = A&/A. The 
classical limit {% — ► 0) of quantized chaotic system (d > 1) is characterized by A <C Aj, 
and hence b 3> 1. But in one dimensional systems id = 1) we do not necessarily have 
this separation of energy scales. For some typical perturbations b — 0(1). We shall 
explain that in such case there is no FGR regime in the theory. A necessary, but 
not sufficient condition for having FGR regime in one-dimensional systems is to have 
6^1. This means that "small features" should characterize the phase space manifolds 
that support the perturbed (or evolving) quantum mechanical states. 

The analysis of one dimensional (d = 1) systems is highly non-universal, but 
typically allows analytical calculations that go well beyond the leading semiclassical 
approximation. In particular we demonstrate the capabilities and the limitations of 
the uniform approximation |12l 

The outline of this paper is as follows: In Section we define some simple 
prototype models. In Section[3]we define the main objects of the studies which are the 
LDOS and the survival probability V(t), In Section^ we discuss the issue of quantum- 
classical correspondence (QCC) and make the distinction between "semiclassical" and 
"perturbative" approximations. In Sections[S]and[S]we discuss the notion of "regimes" 
in the context of LDOS studies. In Sections [7| and [S] we discuss the existence of LDOS 
regimes in ID within the framework of a uniform approximation. In Section [§] we 
discuss the implication of the strong localization effect. In Section ^]] we extend the 
general consideration into the context of survival probability studies. Conclusion and 
final remarks are summarized in Section ITTI The appendices are an integral part of 
the Paper. They contain details of derivations and have not been integrated into the 
main text in order to simplify the presentation of the physical picture. 

X In section 5 we explain that perturbation theory can lead to either no decay, or Gaussian decay, or 
Wigner type decay. The latter can be regarded as the outcome of FGR transitions. 
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2. Definition of the models 

In this paper we are going to consider one-dimensional model systems. The simplest 
is the deformable harmonic oscillator: 

H(Q,P;x) = ±(xP) 2 + \{Q/x) 2 . (2) 
where (Q, P) are the canonical coordinates, and a; is a constant parameter. The energy 
surfaces H(Q, P; x) — E are ellipses in phase space. We shall assume that E 1. 
Without loss of generality we shall regard x = xq = 1 as the "unperturbed" value of 
the parameter x. Later we shall assume that Sx = (x — xq) is small in a classical sense. 
Namely Sx <C 1. Then it is possible to write the Hamiltonian as H = Hq + SH, where 

Ho = H(Q,P;x =l) = i(P 2 + Q 2 ), (3) 



SH = Sx , 
ox 
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= Sx (P 2 - Q 2 ) (4) 

=i 

Note that the perturbation Sx is allowed to be large in a quantum mechanical sense: 
many levels can be "mixed" by the perturbation SH. 

A more complicated case is to consider a particle in a ring (= one dimensional 
box with periodic boundary conditions). The Hamiltonian is of the general form 

H(Q,P;x) = ^-P 2 +x V(Q) (5) 
zm 

Without loss of generality we can take m = 1 as the mass of the particle, and L = 2ir 
as the perimeter (length) of the ring. The parameter x controls the "height" of the 
potential landscape. Naturally we shall regard xq = as the unperturbed value of x. 
We shall consider the case of a single "bump" where 

V(Q) = V exp [- (Q - Q f /2£ 2 \ (6) 

It is implicitly assumed that icL. The Fourier components of this potential are 
iV^fc)! = Vq£ exp [— {k£) 2 /2\. Hence the non-vanishing (\k£\ < 1) Fourier components 
satisfy 1^(^)1 w Vq£. If we have many bumps, then we have a "disordered" potential 

V(Q) = ^(±random)F exp[-(Q-Q Q ) 2 /2^] (7) 

a 

Here we implicitly assume that the bumps are non-overlapping and randomly 
distributed along the ring, such that the correlation function (V(Q + r)V(Q)) is 
characterized by the correlation length £. Consequently the non-vanishing (\k£\ < 1) 
Fourier components of the potential satisfy |t^(fc)| « V Q £x y / L/£. Note that the phases 
of the Fourier components look "random" , unlike the case of a single bump. 

All models that we have introduced are of the generic form H = Ho + SH. The 
representation in the basis that is determined by the unperturbed Hamiltonian is 

H * — > E + SxB (8) 
where E = {E n } is a diagonal matrix consisting of the eigenenergies of the unperturbed 
Hamiltonian. With scaled parameters, such that h = 1, the mean level spacing of the 
eigenenergies is A = 1 for Hamiltonian @ and A = \J2E for Hamiltonian (JSJ. The 
matrix B corresponds to the perturbation. It is a banded matrix. The bandwidth 
is b. It is defined such that 2b is the number of the levels which are coupled by the 
perturbation. For the Hamiltonian J2J) we have b — 1, because only neighboring levels 
of the de-symmetrized Hamiltonian (\n — m\ = 2) are coupled by the perturbation SH. 
For the Hamiltonian © the bandwidth is b = L/£. See Section[7]for details. Thus for 
the latter model we may have b ^> 1 . 
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3. Definitions of P(n\m) and V(t) 

Let \n(x)) and E n (x) be the eigenstates and the corresponding eigenvalues of the 
Hamiltonian H.(Q, P;x). We define the "parametric" kernel 

P(n\m) = \{n(x)\m(x ))\ 2 = trace(p„p m ). (9) 

Note that for Sx — we have P(n\m) = <5„ m . Given a reference state |to(:eq)) this 
kernel can be regarded as a probability distribution with respect to n. The LDOS is 
just a scaled version of this kernel, namely 

P{lu) = ^P(n\m)2w6(cj-[E n (x)-E m (x Q )]} (10) 

n 

One important measure that characterizes the LDOS is its variance: 

6E 2 = $^P(n|m) (E n -E m f (11) 

n 

In the next section we shall explain that the variance has a special role in the theory 
of the LDOS. 

The survival probability of the state \m(xo)} for evolution which is generated by 
the perturbed Hamiltonian H(Q, P; x) can be written as V(t) = |.F(t)| 2 , where 

F(t) = (m(xo)\e- lt ( n - E ^\m(x )) (12) 

In the above definition we have taken E = E m as the natural reference for the energy. 
With this definition we see that F(t) is just the Fourier transform of the LDOS. We 
note that in more complicated scenarios a simple Fourier transform relation between 
F(t) and P(ll>) does not exist UJ. 



4. Quantal-classical correspondence (QCC) 

The classical limit/approximation for the LDOS kernel P(n\m) is obtained by taking 
for the Wigner function p n (Q, P) « p%(Q, P), where p%(Q, P) cx 6 [H(Q, P; x) - E n (x)] 
is the corresponding classical microcanonical distribution. For the deformable har- 
monic oscillator one obtains | |Appendix A| ): 

P*(n\m) » - — 1 (13) 

tt y/4(5x/x)*E* - (E n - E m y 

The dispersion (square root of the variance) is 

SE cl = V2^jE (14) 

What about the quantum mechanical LDOS? The simplest limit that can be 
considered is first order perturbation theory (FOPT). One obtains ( |Appcndix B| ): 

1 fSx x " 



P^(n\m) « S nm + - ^—Ej S {n _ mlt2 [FOPT] (15) 

This result is very different from the classical result. Therefore we say that there is no 
"detailed QCC" . On the other hand one easily calculates the variance. One obtains 
SE = SE cl . We call the latter equality "restricted QCC". 

It is possible to write an exact analytical expression for the quantum mechanical 
LDOS (see |Appcndix A| ). But this expression is not very useful since its complexity 
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makes it virtually impossible to extract the simple limits in various regimes. It is more 
illuminating to obtain a uniform approximation l |Appendix C| ) 

-l 2 



P uniform / i _ \ 
(n\m) 



J(m-n)/2 ( E 



X 



(16) 



In figure^wG demonstrate that the uniform approximation is almost indistinguishable 
from the exact result for any value of Sx. One can verify that this approximation 
reduces to the FOPT result (|15|l in the parametric regime Sx <C (x/E). Disregarding 
an oscillatory component it reduces to the semiclassical result 1|13|) in the parametric 
regime 8x 3> (x/E). The two regimes are separated by the quantum mechanical 
parametric scale Sx plt = (Tiuj osc / E)x. For the convenience of the reader we have 
reverted here to non-scaled units (with our scaling oj osc — 1 and fi = 1.) 

The example above demonstrates some general observations regarding the LDOS. 
On the one hand we have restricted QCC, which means 6E — 6E cl . It can be proved 
that this type of QCC holds in general . The proportionality 8E cx Sx is guaranteed 
by the classical linearization condition, which is a fixed assumption in this paper (in 
the present example it means Sx x). On the other hand we do not have in general 
detailed QCC, which means that the approximation P(n\m) ~ P cl (n\m) holds only 
in a specific parametric regime. In the above example we have only two parametric 
regimes: 

• The "perturbative regime" (Sx <C Sx prt ) 

in which (in this example) FOPT can be used. 

• The "non-perturbative regime" (Sx > 5x pTt ) 

in which (in this example) the semiclassical approximation can be used. 

Hence upon quantization there is only one parametric scale in the theory (Sx plt cx Ti). 
This parametric scale marks a border between the perturbative and the non- 
perturbative regimes. 



5. Regimes in case of chaotic systems 

In case of chaotic systems the generic case is having three regimes. The "perturbative 
regime" is subdivided into a "FOPT regime" and a Wigner/Lorentzian/core-tail/FGR 
regime. In the latter case the various names mean the same, and we shall use from 
now on the term "Wigner regime" . The FOPT regime is defined as the parametric 
regime where we can use first order perturbation theory (FOPT): 

P^(n\m) « S nm + fa 2 |B "" t|2 [FOPT] (17) 

The condition for the validity of this approximation is Sx <C Sx c , where Sx c — A/er. 
Here A is the mean level spacing, and a is defined as the RMS value of the "in-band" 
off-diagonal elements of the B matrix. 

Outside of the "FOPT regime" we can still use first order perturbation theory for 
the tails of the LDOS ^Uj- This leads to a generalized Wigner's Lorentzian (a more 
meaningful name is "core-tail structure" ) : 

P prt (n|m) a Sx 2 !*?nm]_ [Wigner] (18) 

[T(5x)/2] 2 + (E n -E m f 

The width parameter T(Sx) is determined so as to have proper normalization 
(J2 n E(n\rn) = 1). For strongly chaotic systems, for which the band profile is quite 




Figure 1. The local density of states for E=100. (a) For (8x/x)E = 0.5 
(perturbative regime), (b) For (8x/x)E = 2 (intermediate regime), (c) For 
(8x/x)E = 20 (scmiclassical regime). In each panel the thick dashed line is 
the exact result, the solid line is the uniform result, the dashed-dotted line is the 
classical result, and the solid circles are the perturbative result. The classical 
result is divided by 2 for a reason which is explained at the end of|Appcndix A| 
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flat, it follows that the "core" width is T w (Sx/Sx c ) 2 A. The "core-tail" approximation 
makes sense as long as we have separation of energy scales T(Sx) <C A&, where A& = 6A 
is the bandwidth in units of energy. This translates into the condition 8x <C Sx pTt 
where 

Sx plt = Vb Sx c . (19) 

The "core-tail" approximation (|18fl can be regarded as the outcome of 
perturbative summation of diagrams to infinite order: The FOPT diagrams are re- 
iterated, while the interference terms are neglected. In time domain analysis this 
corresponds to a Markovian approximation for the survival probability V(t), leading 
to an exponential decay. Hence we say that there is a perturbative regime that includes 
the FOPT sub-regime, and the Wigner sub-regime. In the Wigner sub-regime we need 
all orders of perturbation theory leading to FGR transitions and Wigner decay. On 
the other hand in the FOPT sub-regime there is no decay. This can be easily deduced 
by Fourier transforming the LDOS which is associated with Ea. p7|l : Up to negligible 
first order correction, the survival amplitude F(t) comes out as a pure phase factor, 
whose absolute value squared is V(t) 1. It is appropriate at this point to make 
a connection with recent fidelity studies 4 . There, it is customary to consider the 
dynamics of a wavepacket that is a superposition of many eigenstates. Consequently, 
in fidelity studies, there is an effective averaging over the survival amplitude F(t), 
leading (in the FOPT regime) to a slow Gaussian decay which is trivially related to 
the statistics of the first order correction to the eigenenergies. 

We can summarize this section by stating that in case of a generic quantized 
chaotic system we have three regimes as follows: 

• The FOPT regime (Sx < 5x c ) 
where Ea. l(T7|) can be trusted. 

• The Wigner (perturbative) regime (5x c <C Sx <C Sx VIt ), 
where Ea. l(T%| can be trusted. 

• The non-perturbative regime (Sx > <5£ prt ), 

in which (typically) the semiclassical approximation can be used. 

Hence upon quantization there are two parametric scales in the theory, which are Sx c 
and Sx plt cx h. These parametric scale mark the borders between the regimes. 

6. Regimes in case of ID systems 

By Weyl law we know that the mean level spacing for d-dimensional system is A oc Ti d . 
On the other hand the bandwidth A& cx Ti is inversely proportional to the period, or to 
the correlation time of the dynamics. For chaotic systems (with d > 1) we generally 
have b — A&/A ^ 1. But for d = 1 systems we can have b = 0(1), as in the 
example of equation @. Therefore in the latter case Sx pTt ~ 5x c , and we do not 
have a "Wigner regime" . This means that the LDOS cannot have a Lorentzian-likc 
structure, and consequently the survival probability V(t) = \F(t)\ 2 cannot have an 
exponential decay. While this latter statement strictly holds for the specific scenario 
which has been defined by equation (fl^|l . it typically holds also in more complicated 
circumstances ( "fidelity" studies) . Thus we conclude that the FGR picture cannot be 
applicable to the analysis of the dynamics unless 6 3> 1. 

At this stage of the discussion we can say that in order to have three regimes 
(FOPT, Wigner, semiclassical) in the theory of one-dimensional systems, we have to 
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consider models where the perturbation is characterized by a large bandwidth (6 ^> 1). 
In the next sections we would like to further discuss the following: 

• A uniform approximation can be applied in order to address both the perturbative 
and the non-perturbative regimes. 

• We can have b 3> 1 by considering potentials that have small (or sharp) features. 
(For example "bump" or "disorder"). 

• It is possible to get a Lorentzian-like LDOS from the uniform approximation, but 
6 3> 1 is not a sufficient condition. 

• The uniform approximation does not take into account the effect of strong 
localization. 

• Three LDOS regimes can be observed in case of disordered potential in spite of 
the strong localization effect. 



7. Perturbations with large bandwidth 

For any one-dimensional system we can find the action-angle variables, in which the 
problem becomes essentially equivalent to the problem @| of a particle in a ring. 
Consider for example the deformable harmonic oscillator @: We already saw that 
all the essential physics of the perturbed eigenstates can be obtained via a uniform 
approximation which is based on action-angle variables description of the system. The 
Hamiltonian l|C.3j) of a deformable harmonic oscillator in action-angle variables, and 
the Hamiltonian Q of a particle in a ring, are similar as far as the semiclassical 
treatment is concerned. Having a quadratic rather than a linear dispersion relation is 
not an essential difference. 

More generally, we may encounter circumstances where the perturbation creates 
small phase space structures. This is the new ingredient (compared with the 
deformable harmonic oscillator) that we are going to consider in the following sections. 
The problem J^J of a particle in a ring constitutes a prototype example of having such 
small phase space structures. The perturbation V(Q) is characterized by a different, 
much smaller scale (£), compared with the size of the system (L). This means that 
we have large bandwidth (b 3> 1) rather than b = 0(1). 

In the absence of a perturbation (a? = Xq = 0), the eigenenergies of a particle 
in a ring are E n = p^/2m, where n is an integer index, and p n — (2Trh/L)n. The 
perturbation matrix is related to the Fourier components of the potential: 

B„ m = jV{j) n - Pm) (20) 

The expressions for V(k) in case of either the "bump" © or the "disorder" (JJJ), were 
given in Section [5] Substitution in Q20|l allows to determine the bandwidth: 

A = 2Trhv E /L (21) 
A b = 2Trhv E /l (22) 
b = Lji (23) 

where ve = \/2E/m. One difference between the two models is related to the coupling 
parameter <j, which has been defined in the beginning of Section [5] as the RMS value 
of the in-band off-diagonal elements: 

a = (£/L)V for the bump (24) 
a = {£/L) 1/2 V Q for the disorder (25) 
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This leads to 

Sx bum P = A/ff = {Ab/VQ) (26) 
^disorder = A/ff = ( Ai /y )/V6 (27) 

However, the difference in cx, and hence in Sx c is not the significant difference between 
the "bump" potential and the "disorder" potential. The significant difference is related 
to the statistical properties of the B matrix, tn the case of "disorder" the matrix 
elements look "random" , which is not the case for a single bump. 

In both cases, of having either single bump or disorder, the FOPT regime is 
Sx <C Sx c - What do we have beyond FOPT? More specifically, the question is whether 
we have, as in the theory of chaotic systems, a distinct parametric scale 5x pr t = VbSx c 
that distinguishes between a "Wigner regime" and a "semiclassical regime". In the 
next section we shall try to address this question within the framework of the uniform 
approximation. 



8. Regimes within the framework of the uniform approximation 

There is a general semiclassical procedure that associates wavefunctions of integrable 
systems with phase space manifolds. The traditional implementation of this procedure, 
in case of one dimensional systems, is known as "the WKB method" . The 
WKB method is problematic near turning points. This problem can be solved 
by "uniformization" of the solution. The simplest point of view regarding this 
"uniformization" is obtained by using action-angle variables. This leads to a 
Hamiltonian that looks like that of a particle in a ring. In general the dispersion 
relation can be different (not quadratic as in ©), but we shall see that this is not an 
important difference for the issues under study. 

The eigenstates of Hamiltonian JSJ are supported by the manifolds H(Q, P; x) = 
E n . An alternate (explicit) way to describe a given manifold is P — p n {Q', x), where 



Pn{Q\x) = v'Zm [E n -SxV{Q)} r» p n - 5x ft k(Q) (28) 

and k(Q) = V(Q)/Tlve- In the following we use units such that ft — 1. The 
semiclassical formula for the corresponding eigenfunction is 

(Q\n(x)} = ^ ex P^ j\n(Q';x)dQ'^J (29) 

This can be regarded as a special case of <|C.5|) . Above we approximate the classical 
pre-exponential prefactor by 1/ \fh. This is legitimate because a fixed assumption of 
this Paper is that we are considering high lying eigenstates, and assume classically 
small perturbations §. We already saw, in the context of the deformable harmonic 
oscillator (see the end of | Appendix C| , that the numerical error which is associated 
with this approximation is insignificant. The essential physics of having various 
"regimes" is not related to this approximation. 

The overlap of the semiclassical wavefunctions is given by the integral 

{n(x)\m(x )) = -jr J dQ exp -i ( (p n - p m )Q - 6x J k{Q')dQ^j (30) 

In the following subsections we discuss the consequences of this expression in case of 
either bump or disorder. 

§ It can also be regarded as an example of a more general semiclassical perturbation approximation 
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8.1. bump case 

In case of the bump, k(Q) has an amplitude k$ 
total phase variation in (|3UII is 



Vq/Uve over a spatial scale t. The 



^bump = Sxx k £ (31) 

This phase variation is in fact the phase space area of the bump. So we have two 
possibilities: Either 50bump <C 2ir or <$</>bump 2ir. This can be shown to be equivalent 
to either 8x <C 8x^ ump or 8x ^> fe]? ump respectively. In the former case it is easy 
to recover the FOPT result l|17l) . One should simply put the perturbation off the 
exponent, and then do integration by parts. For n ^ m it leads to 

2 



P(n\m) 



T / dQ- 



-i(Pn-Pm)Q 



Pn ~Pv 



Sxk{Q) 



Sx 2 



B, 



E n — E 

I) 



(32) 



The other possibility (54>bump 3> 2ir) guarantees the validity of the standard 
semiclassical approximation. In such case we cannot put the perturbation off the 
exponent, but instead we can make a stationary phase approximation. We shall not 
dwell further on details because it is a standard textbook procedure. 



8.2. disorder case 

We see that for a simple bump we have either FOPT or semiclassical approximation. 
So we have just two regimes, in spite of the fact that b ^> 1. This means that 
large bandwidth is not a sufficient condition for having three parametric regimes. So 
let us try to make things more complicated by considering a disordered potential 
with many bumps. Equation (|31|l still describes the phase variation over a single 
bump. This means that 50bump ^S> 2ir is still the relevant condition for a semiclassical 
approximation! What about first order perturbation theory? The total phase variation 
in (|30|l for many bumps is 



5x x Vb x k £ 



(33) 



(note that b is essentially the number of bumps involved). Consequently the validity 
condition for first order perturbation theory is 50disordcr <C 2ir, which can be easily 
converted into Sx <C fe^ lsordcr . 

The considerations of the previous paragraph imply that for disordered potential 
we have three regimes. In the intermediate regime Sx c < fe < 5x VTt , we have 
^disorder > 2lr while 6$ nmp < 2tt. Consequently we can use neither FOPT, nor 
the semiclassical approximation. This is the Wigner regime where we expect to find a 
Lorentzian-like LDOS. Let us demonstrate that indeed a Lorentzian-like LDOS can be 
obtained from the uniform approximation (|30|l . For this purpose we average P{n\m) 
over realizations of the disorder: 

rQ 2 



P{n\m) 



1 

1 

1? 



dQidQ2 e 



-i(p„-p m )(Q 2 -Qi) 



exp 



idx / k{Q')dQ' 



rdcr 



dQidQ 2 e 



dr e- l(p ^ p ^ r 



exp 



1 f—Y 

2 \hv~E~J 



Q2 i-Q 



{V{Q')V{Q"))dQ'dQ" 



exp 



( Tiv E ) 



(34) 
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The last integral is the Fourier transform of an exponential, leading to the Lorentzian 
LDOS as defined in ijlgjl. 



9. Strong localization effect 

Still we have to address the question whether we can trust the uniform approximation, 
which is based on the WKB wavefunction (j29(l . The answer is known to be negative 
in case of disordered potential. The WKB approximation does not take into account 
backscattering, which is responsible for the strong localization effect in ID disordered 
potential. It is well known ^3 that the localization length is equal (in ID) to twice 
the mean free path. Up to a numerical prefactor the Born approximation estimate is 

1 la 2 „ 2 f 5x\ 2 I f Sx \ 2 1 



L loc - hv E A X5X ~ [sxj L {sx p J I (35) 

The condition for not being affected by the strong localization effect is L loc 3> L, which 
leads to Sx <C Sx c . Thus it follows that only in the FOPT regime we can ignore the 
strong localization effect. The strong localization effect cannot be ignored neither in 
the Wigner regime nor in the semiclassical regime. 

The above Born approximation assumes L lac 3> I. This condition breaks down 
if Sx > Sx prt . Recall that we also assume that Sx is small in the classical sense 
(SxVo <C E). The two inequalities are consistent if and only if the de Broglic 
wavelength of the particle is much smaller compared with I. In this regime we can 
analyze the localization using the well known transfer matrix approach: Each bump 
has some transfer matrix, and the random distance between the bumps provides the 
phase randomization which is assumed in "combining" adjacent transfer matrices. 
Denoting the average transmission of a "bump" by g, one gets 

J- » Hl/g) x j (36) 

Thus even in the Sx > Sx VIt regime we can have a very long localization length 
(Li oc 3> £), which is in fact consistent with the naive expectation. 

The implications of the above discussion are, that in spite of the strong 
localization effect, it is still meaningful to distinguish between three parametric regimes 
(FOPT, Wigner, semiclassical). We just have to remember that the wavefunction is 
not ergodic in real space, so the role of L is taken by L loc . 



10. The survival probability V(t) 

We turn to discuss the calculation of the survival probability QJ for the specific 
"wavepacket dynamics" scenario that has been defined in Section We have the 
following five strategies of calculation: 

• Uniform approximation (which is essentially exact) 

• Time domain classical approximation 

• Energy domain classical approximation (~~» LDOS ~> Fourier transform) 

• Time domain perturbation theory 

• Energy domain perturbation theory (~» LDOS ~> Fourier transform) 

It can be shown that in typical circumstances the two versions of perturbation 
theory give in leading order consistent results. This means that we can write a 
perturbative (essentially first order) result that can be trusted for sufficiently short 
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times (for any perturbation 5x) or for sufficiently weak perturbation (for any time t) . 
As an example we consider the deformable harmonic oscillator. Taking equation (|15fl 
with appropriate second order compensation of normalization, we get after Fourier 
transform, 

V prt (t) » 1 ~ 2 (^) E 2 sm 2 t (37) 

In a strict time domain FOPT we get only t 2 time dependence, while in a strict energy 
domain FOPT we do not get the correct normalization (P prt (0) = 1). Still we are able 
to get one consistent result. The situation is different with the classical approximation. 
Here time domain and energy domain calculations do not give the same result. As 
an example we consider again the deformable harmonic oscillator. The calculation of 
overlap between pf and Pq is simpler in action-angle variables, but otherwise it is 
similar to the calculation in | Appendix A| leading to 

-l 



V c \t) 



Sx 

2ir — Esmt 
x 



(38) 



The energy domain classical approximation is obtained by squaring the Fourier 
transform of equation (fTT)fl using the classical approximation lfO| of P cl (n\m), leading 
to 



■ycl,E 



(*) = 



, 8x 
Jo ( 2—Et 

x 



(39) 



Although there is no simple exact expression for V(t), there again exists a very simple 
uniform approximation which is remarkably accurate in all regimes and therefore we 
can regard it as "exact," namely 

2 



3 uniform 



-/,: j 2— Esint 



(40) 



It is derived in |Appendix D| by using semiclassical expressions for the initial and 
evolved states, while calculating the overlap exactly rather than by the stationary 
phase approximation. 

Results <|57J) - (EiU|) are graphically displayed in figure[21 Approximations 1(3*7 )1 -1(55 )1 
can be regarded as various limits of the uniform approximation 1)40)1. It is easily seen 
that (fit])) reduces to the perturbative result ||S7J whenever 2(Sx/x)Esin(t) <§; 1. So 
as expected the perturbative result can be trusted for either small time t or for small 
perturbation Sx/x. It is also easy to see that the uniform result 1)40(1 reduces to the 
energy domain classical result l|39|) for t <C 1. The relation of 140)1 to the time domain 
classical result 1)38)1 is more subtle: For large perturbation the two expressions agree 
in an asymptotic sense, and either time smoothing or energy averaging is required in 
order to demonstrate this agreement (see figure Et). 

An important message of this section is that the discussion of regimes in the 
context of V(t) is in one to one correspondence with the discussion of LDOS regimes. 

11. Conclusions and final remarks 

In general, for more complicated "wavepacket dynamics" scenarios, which may involve 
time dependent Hamiltonian, an explicit reduction of the survival probability problem 
to LDOS study is not possible [SJ. The simplest way to make the Hamiltonian "time 
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Figure 2. The survival probability, using the same parameters as in figure^ (a) 
Perturbative regime, (b) Intermediate regime, (c) Semiclassical regime. In each 
panel the solid line is the uniform result 1401 . the dashed line is the perturbative 
result 1371 . the dashcd-dotted line is the time domain classical result 1381 . and the 
dotted line is the energy domain classical result 1391 . Panel (d) gives a short-time 
detail of panel (c). 
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dependent" is by changing it either once (as in "fidelity" or "Loschmidt echo" studies) 
or repeatedly (as in kicked systems) . The latter possibility opens the way for "chaos" 
in the classical dynamics of one-dimensional (ID) systems. A prototype system for ID 
chaos studies is the kicked rotator: This is a particle in a ring, which is periodically 
kicked by a cosine potential. As in the standard formulation of "wavepacket dynamics" 
it is common to assume an initial preparation which is supported by a manifold 
P = Pm ■ As a result of the kicks the initial manifold becomes extremely convoluted. If 
we want to calculate V(t) we have to trace the evolving manifold with the initial one. 
The calculation of this overlap can be carried out using the uniform approximation 
that we have discussed previously [T^HTHl I?] . 

There is one major difference that makes the "chaotic" scenario that we have 
described above different from the LDOS calculation: The issue of strong localization, 
which has been discussed in Section is no longer relevant. The evolving state is 
"distributed" over the whole manifold ||. 

In view of the above, the calculations that we have presented in this paper, using 
the uniform approximation, are in fact generic. In particular we have explained how a 
Fermi-golden-rule behavior arises within this framework. A refinement of the uniform 
approximation using the "replacement manifolds" approach has been introduced in 
references |12l 113) . 

The reader can be tempted to reach the conclusion that the discussion of 
"regimes" within the framework of the "uniform approximation" can be trivially 
generalized from d = 1 chaotic systems to d > 1 chaotic systems. This is in fact not 
quite correct. In order to explain the difficulty let us consider again LDOS calculation. 
In one-dimensional systems what we have to do is to calculate the overlap of two one- 
dimensional manifolds. We can do this calculation semiclassically if we can trust the 
stationary phase approximation. This leads to the " Ti area" condition: the stationary 
phase approximation is accurate if the phase space area delimited by the two manifolds 
between two stationary phase points is larger than h. What is the generalization of 
this condition in the d > 1 case? Now the manifolds are "surfaces" and they intersect 
along "lines" . So the concept of "stationary points" becomes inapplicable, and also the 
u h area" condition becomes meaningless. The way out of this difficulty 10 is to use 
the Wigner function point of view. Then one realizes that the proper way to formulate 
the semiclassical condition is to say that the separation between the surfaces should 
be much larger compared with the "thickness" of Wigner function. This "thickness" 
is just the "bandwidth" that we have discussed in this Paper. 

Conventionally, semiclassical methods are applied in the "semiclassical" regime 
(large enough perturbation). In the perturbative regime people use perturbation 
theory. A major motivation for the present line of study is to extend the applicability 
of semiclassical methods into the perturbative regime. Most of the calculations that we 
have presented under the heading "uniform approximation" can be reformulated using 
the Wigner function language. This opens the way towards a unified semiclassical 
understanding of "regimes" in case of d > 1 systems. 

|| The only reservation for this statement is the possibility of witnessing a "dynamical localization 
effect" after an extremely long time (known as the "breaktime" ) . For small Ti the existence of a 
"breaktime" will have a negligible effect on the behavior of V{t). In fact we should remember that 
also in non-kicked systems we have a "breaktime", which is just the Hciscnbcrg time. We also note 
that this type of localization is apparently captured by semiclassical methods 
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Appendix A. Calculation of P cl (n\m) and P oxact (n|m) 

We shall demonstrate in this appendix the calculation of P cl (n\m) for the deformable 
harmonic oscillator: 

P cl (n\m) = [^^5[H(Q,P;x 2 )-E n ]5[H(Q,P;x 1 )-E m ] 

dQ ^ 6[H(Q,P;x)-E n ] 

2?r l± X\\P\ 

P=±y / 2E m -(Q/x 1 )2/x 1 

1 f d Q- WW 



where 



J y/2x\E m - Q 2 

l^2(HE m -Qiy 1/2 \f(Q±)r 

1 X\X 2 

n y/{x\E n ~ x\E m ) {x\E m ~ x 2 E n ) 



4 _ 4 / \ 2 



(A.I) 



/(Q) ee ^-^Q 2 + P ) E m - E n . (A.2) 



4 2 

and Q± are the roots of the equation f(Q) = 



r\ _i_ 2(x 2 E m X 2 E n ) /, „> 

Q± = ±XiX 2 \ 4 _ 4 (A.3) 

y x 2 Xi 

for which 

f\Q±) = ±x^ 3 x 2 1 ^2(x 2 E m - x\E n )(x\ - xj) (A.4) 

Equation (|13|l is a simplified version of this result, assuming that Sx <C x. 

It is also possible to obtain an exact result in the quantum mechanical case. The 
explicit expression for the eigenfunction using Hermite polynomials is 

(Q\n(x)) = {TTx 2 )- 1 ' i {2 n 7i\)' 1 ' 2 H n (Q/x)e- {Q/x)2/2 (A.5) 

This leads to 

(n(x 2 )\m(x 1 )) = {■Kx 1 x 2 )- 1 ' 2 (2 n+m n\m^- 1 ' 2 

djQH n {Q/x 2 )H m {Q/x l )e-\ { - x " +x ^ Q2 (A.6) 

Upon squaring one obtains P(n\m). The integral in 1JA.6J1 becomes highly oscillatory 
for high-lying eigenstates in which we are interested, and numerical calculation is 
tricky. To our surprise, this intimidating integral can be evaluated analytically even 
for x\ 7^ x 2 , resulting in a finite sum of terms, which is not very elegant, but very simple 
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to evaluate on a computer. For the exact benchmark in the numerical results presented 
in this paper, we therefore used this analytical expression instead of numerically 
evaluating the integral (|A.6(1 . 

Note that both classically and quantum mechanically the overlap P(n\m) depends 
only on the ratio xij X \. In the classical case n is a real index (E n can have any real 
positive value). In the quantum mechanical case n is an integer index (E n = n + |). 
Due to the reflection symmetry of the Hamiltonian there is an overlap only between 
states with the same parity. The overlap for \n — m\ — odd vanishes. Whenever we say 
that P(n\m) P cl (n\m), it should be interpreted in a coarse grain sense. Therefore, 
for sake of graphical presentation we have plotted P(n\m) versus P cl (n\m)/2. One 
may say that the plotted P(n\m) and P cl (n\m)/2 correspond to a de-symmetrized 
oscillator. 



Appendix B. Calculation of P prt (n|m) 

We shall demonstrate in this appendix the calculation of P prt (n\m) for the deformable 
harmonic oscillator. From first order perturbation theory we know that for m =/= n 

( m ( Xo )\n( Xo +5x))^ { ^ m ^ (B.l) 

Using 



<m|Q 2 |n), (m\P 2 \n) = (± + n)S mn ± - |y (ra + l)(n + 2)5 m , n+2 + ^n{n - l)£ m ,„_ 2 
one obtains 

P^(n\m) « S n , m + \ (e 2 - V\ 5 ln _ ml>2 (B.2) 

where E = (E n + E m )/2. In the text we have presented a simplified version that 
assumes E 3> 1 . 



Appendix C. Calculation of F uniform (n|m) 

The following canonical transformation is used in order to transform the Hamiltonian 
of the harmonic oscillator to action-angle variables: 

Q = \/27cos(» (C.l) 

P = V2Isw((t>) (C.2) 

In the vicinity of xq = 1 it leads to 

H{(j>,I;x) = I [(l/x)cos 2 + xsm 2 (f>] 

w 1-5x1 cos(20)) + 0{5x 2 ) (C.3) 

The canonical transformation from (4>,I) to the action angle variables (</)', I') of the 
perturbed Hamiltonian is derived from a generating function S(<fi,I'). The manifold 
/' = const is determined from the equation 7i (</>, /; x) = const, leading to the relation 
I = [1 + <5x cos(20)] The generating function should satisfy I = dS/d<p. Therefore 
one deduces that 



S{(j),I') = [4> + Sx sin(2<^>)] I' 



(C.4) 
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The semiclassical expression for the wavefunction is 

The nth semiclassical eigenstate corresponds to the substitution I' = E n = n + 1 . For 
technical simplicity it is convenient to calculate the overlap between two perturbed 
wavefunctions [x = x ± Sx/2) . 

{n(x + Sx/2)\m(x -Sx/2)) = / #(0|/"=£ n )*(0|/'=£ m ) 

Jo 

I /-27T 

— — d(f>\/ 1 — (5a; 2 cos 2 20 x exp [i(i? m — E n )(f> — iSxE sin 20] 
2tt Jo 

~ J(,n~n)/2(SxE) (C.6) 

Above 7' and I" correspond to the perturbations ±<5x/2, and £ = (E n +E m )/2. In the 
main text we have reverted to a more general version of this expression that does not 
assume xo = 1. It is important to realize that because of the assumption Sx -C £ the 
pre-exponential factor (which is in fact a "classical" factor) can be neglected in leading 
order. Figure^confirms the remarkable agreement of the uniform approximation with 
the exact result. 



Appendix D. Calculation of V 



uniform 



it) 



For the purpose of calculation it is more convenient to regard the prepared state as a 
"perturbed state" and the evolution Hamiltonian as the "unperturbed Hamiltonian" . 
This is of course equivalent to the presentation in the text upon the replacement 
Sx i > —Sx. 

The initial state \m(xo + Sx)} is characterized by the action I' = E rn = E. It is 
represented as in | Appendix C| by 

(<p\I'=E) ps — exp[i(0 + Sx sin 20)£] 



2tt 

The evolving state is represented by 
-an 



e- %tUa \I'=E) = (<f>-t\I'=E) 



The overlap between the evolving and the initial state is 

r2w 



1 

2^ 



>exp 



■ ( - t + Sx [sin 2(0 -t)- sin 20] J E 



For F(t) as defined in (|12fl we obtain (after the required replacement Sa 

r2n 



F(t) 



1 

2^ 



;exp [iSxEsmt cos(20 — t)] — Jo(2SxEs'mt) 



The survival probability is obtained by squaring this result. 



(D.l) 
(D.2) 



-Sx): 
(D.3) 



References 

[1] Cohen D and Heller E J 2000 Phys. Rev. Lett. 84 2841 

[2] Heller E J 1991 Chaos and Quantum Systems, ed M-J Giannoni et al (Amsterdam: Elsevier) 
[3] Kallush S, Segev B, Sergeev A V and Heller E J 2002 J. Phys. Chem. A 106 6006 
Sergeev A V and Segev B 2002 J. Phys. A: Math. Gen. 35 1769 



Survival probability and LDOS for ID systems 



18 



[4] Jalabert R A and Pastawski H M 2001 Phys. Rev. Lett. 86 2490 

Cucchietti F M, Pastawski H M, Jalabert R 2000 Physica A 283 285 
J.V. Emerson, Phys. Rev. Lett. 89, 284102 (2002). 

Cucchietti F M, Pastawski H M and Wisniacki D A 2002 Phys. Rev. E 65 045206 
Cucchietti F M, Lewenkopf H, Mucciolo E R, Pastawski H M and Vallejos R O 2002 Phys. Rev. 
E 65 046209 

Benenti G and Casati G 2002 Phys. Rev. E 65 066205 
Cerruti N R and Tomso vic S Phys. Rev. Le tt. 88 054103 
Prosen T 2001 Preprmt |quant-ph/0106149| 
Prosen T and Znidaric M 2001 J. Phys. A 34 L681 

Eckhardt B 2003 J. Phys. A 36 371 

Silvestrov P G, Tworzydlo J and Beenakker C W J 2002 Preprmt |nTin" CD/0207002 

[5] Jacquod P, Silvestrov P G and Beenakker C W J 2001 Phys. Rev. E 64 055203 

[6] Wisniacki D A and Cohen D 2002 Phys. Rev. E 66 046209 

[7] Vam'cek J and Heller E J 2003 Preprint quant-ph/0302192 

[8] For a pedagogical presentation, including references, see references |5] and 1101 . which can be 

downloaded from http://www.bgu.ac.il/~dcohen 
[9] Cohen D 2002 "Driven chaotic mesoscopic systems, dissipation and decoherence" in Dynamics 

of Dissipation: Proceedings of the 38th Karpacz Winter School of Theoretical Physics ed P 

Garbaczewski and R Olkiewicz (Springer) 
[10] Cohen D 2000 "Chaos, dissipation and quantal Brownian motion" in New Directions in Quantum 

Chaos: Proceedings of the International School of Physics Enrico Fermi Course CXLIII ed 

G Casati, I Guarneri and U Smilansky (Amsterdam: IOS Press) 
[11] Cohen D 1999 Phys. Rev. Lett. 82 4951 
[12] Vam'cek J and Heller E J 2001 Phys. Rev. E 64 026215 
[13] Vam'cek J and Heller E J 2003 Phys. Rev. E 67 016211 
[14] Cohen D and Kottos T 2001 Phys. Rev. E 63 36203 
[15] Miller W H 1974 Adv. Chem. Phys. 25 69 
[16] Miller W H and Smith F T 1978 Phys. Rev. A 17 939 
[17] Richter K, Ullmo D and Jalabert R A 1996 Phys. Rev. B 54 R5219 
[18] L. Kaplan, Phys. Rev. Lett. 81, 3371 (1998). 



